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Abstract 

A method for classifying n-species reaction-diffusion models, admitting 
shock solutions is presented. The most general one-dimensional two- 
species reaction-diffusion model with nearest neighbor interactions admit- 
ting uniform product measures as the stationary states is studied. Satisfy- 
ing more constraints, these models may experience single-shock solutions. 
These models are generalized to multi-species models. The two-species 
models are studied in detail. Dynamical phase transitions of such models 
are also investigated. 
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1 Introduction 



The stochastic modehng of systems is a useful method for studying the problems 
in non-equilibrium statistical physics. Reaction-diffusion systems are stochas- 
tic models which can be used to study the evolution of interacting particle 
systems. Extensive researches have been done on one-dimensional reaction- 
diffusion systems, some of which belong to the emergence and evolution of 
shocks whose positions perform random walks, i.e. density discontinuities which 
are randomly on the move. The simplicity of the asymmetric simple exclusion 
process, including just diffusive processes, provided researchers with a suitable 
ground to take first steps in exploring the dynamics of shocks [J-E] . In analogy 
to ASEP some interesting models have been introduced, for example a driven 
diffusive two-channel system [B] or bricklayers' model, which is a model without 
exclusion but yet uncorrelated Ol] . To take into account the systems including 
interacting processes, for instance, shock formation in driven diffusive systems, 
containing homogeneous creation and annihilation of particles has been investi- 
gated [QlfTO] and recently more complicated systems have been described [I2l - [l6] . 

It is known that the ordinary Glauber model on a one-dimensional lattice 
with boundaries at any temperature, shows a dynamical phase transition |18j . 
The dynamical phase transition is controlled by the rate of spin flip at the 
boundaries, and is a discontinuous change of the derivative of the relaxation time 
towards the stationary configuration. In [19' , using a transfer matrix method, it 
is shown that a one-dimensional kinetic Ising model with nonuniform coupling 
constants may exhibits a dynamical phase transition. Other phase transitions 
induced by boundary conditions have also been studied ( [20H22] for example). 

We are interested in the works aim to introduce new solvable models. In [TT| 
a single-species model with nearest neighbor interactions on a one dimensional 
lattice with open boundaries has been considered. It was shown that there are 
three families of models with traveling wave solutions; the asymmetric simple 
exclusion process (ASEP), the branching-coalescing random walk (BCRW) and 
the asymmetric Kawasaki-Glauber process (AKGP). A classification of single- 
species models with three-site interactions and special choice of symmetries has 
been studied in [12] . Recently some efforts have been made to obtain the models 
on a lattice with two types of particles and nearest neighbor interactions. In |13) 
a model with diffusion and exchange processes has been studied. In [14 a model 
with a degenerate conservation law and PT invariance (invariant under the ap- 
plication of time reversal and space reflection) has been discussed. Another 
class of two-species models with a non-degenerate conservation law has been 
presented in [15_. Besides the straightforward calculation of the master equa- 
tion there is an alternative approach in dealing with reaction-diffusion problems 
which is the so-called matrix product formalism. It is an algebraic method that 
takes advantage of non-commutative operators instead of probabilities. In jl6) 
it is assumed that the density of A particles in a given site is proportional to 
the density of B particles at the same site. Using matrix product formalism, 
they have found three three-states models which are basically two-state systems. 
One of these models is the generalization of [14] . 
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In this article there is an attempt to present a method for classifying n- 
species particle systems, admitting single shock solutions. In section 2, after 
a brief review of formalism, a two-species reaction-diffusion model on a one- 
dimensional lattice with boundaries is introduced. In section 3, the most gen- 
eral two-species models with nearest neighbor interactions, admitting uniform 
product measures as the stationary states are studied. For such models reaction 
rates should satisfy some constraints. These conditions are obtained. In section 
4, a single-shock measure is introduced and a classification method is presented. 
Then, it is generalized to multi-species models. In section 5, dynamical phase 
transitions of some two-species models are investigated. Finally, a discussion on 
the results of this article and those of previous ones is presented. 



2 Formulation 

Two-species reaction-diffusion models with nearest-neighbor interactions on a 
one-dimensional lattice with L sites are studied. There are two types of particles 
denoted by A and B. The processes arc exclusive which means each site is 
either empty (will be shown by 0) or occupied by at most one particle, A or 
B. The empty state is denoted by |0), and the occupied state by particle A 
(B) is represented by |A) (|B)). The 3-dimensional vector space for each site is 
spanned by 



10) = , |A) = 1 , |B) = . (1) 






The 3''"-dimensional vector space of the lattice is given by the tensor product of 
the single-site vector spaces 



L 

The state vector of the system is 

3^ 

\p)^ = Y,PM\n), (3) 

r,= l 

where |?7) is the basis vector of the lattice and P{r],t) is the probability of 
finding the system in state r] in time t. Continuous-time Markovian evolution 
is the master equation 

i^^Hip). (4) 

The non-diagonal elements of the generator H are the transition rates Ljji, the 
transition rate from state i to j, and its diagonal elements are the negative sum 
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of the non-diagonal elements of their own columns. So the sum of each column 
is zero, and consequently 

{S\H = 0, (5) 
where {S\ is a row vector whose all elements are equal to one. Let us denote 

For a two-species model the local hamiltonian hk,k+i is a 9 x 9 matrix, 
which acts on sites k and k + 1, and contains 72 two-site transition rates. If the 

singlc-sitc states arc respectively assumed empty (0), occupied by A and B, the 
two-site states respectively become: 

1) 00 4) A0 7) B0 

2) 0A 5) AA 8) BA 

3) 0B 6) AB 9) BB (6) 

Let us assume open boundary conditions, where the particles can enter and 
leave the lattice from the first and last sites with the following rates 

Qi2 : A ^ 0, a2i : A 

ai3 : B 0, aai : B (7) 

and change to each other with the rates 

a23 : B ^- A, : A ^- B (8) 

The hamiltonians Bi and Bl are for the left and right boundaries and are in 
the form 



-{0t21 +"31 ) 



-Bl/L = 



21 "T "31 ) "l2 "l3 



"21 -("12 +"32 ) "23 ' 



23 



where the indices £ and r stand for left and right. So the hamiltonian for a 
one-dimensional L-site lattice with nearest-neighbor interactions and single-site 
interactions at the boundaries is in the following form 



H = Bi^ + J2 l®^'^-^) hk,k+i O l®^^-'^-^) + 1^(^-1) (g) Bl (10) 



k=l 



1 denotes a 3 x 3 unit matrix. We assume the interactions to be homogenous, 
so hk,k+i is the same for all sites. We also assume the interactions to be time 
independent. 



3 Uniform product measures as stationary states 

Let us consider the state vector of the system to be uncorrelated. Then the 
state of the system is the tensor product of the single-site state vectors 

\P) = \Pi)^\P2)®---®\Pl) (11) 
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where \Pk) is the state vector of site k. If the occupation probabihty of particle 
A in site k is Uk and the occupation probabihty of particle B in site k is bk, the 
probability of being empty becomes (1 — (ofe + 6fc)). Consequently 



Pk) = flfe (12) 




where not only < Ofc , bk < I but also < + 6^ < 1. We also assume 
that the densities are uniform, and the state vector of each site is shown by \u). 
Then the state vector of the lattice becomes 

\p) = (13) 

which is a uniform product measure. The state vector \P) is the stationary state 
of the system if 

H\P) = 0, (14) 

which gives 

0, (15) 

One can expand h\u) ® \u) in the following form 

2 2 

h\u)(i)\u) = A\u)(E)\u)+ Bf,^\x^)^\x^) + J2{(^t^W)®\xt^)+Di,\Xf,)®\u)}, 

/i,!' — 1 fl—1 

(16) 

where A, B^^i,, C^, and are constant, and the two vectors \xi) and 1x2) 
together with |m) form a linearly independent set. Hence, substituting the above 
expansion into the equation (jlSp leads to a combination of linearly independent 
terms. Thus, one deduces 

A= B^, = 0, C^ + D^=0. (17) 

Therefore, for an infinite lattice or a periodic one, where there is no boundary 
term, equation (|16l) recasts to 

h\u)(S)\u) ^ \u) (S) \x) - \x) \u) , (18) 

where the vector |a;) is a linear combination of \xi) and \x2)- This gives some 
constraints on the reaction rates. For a lattice with boundaries, \u)'^^ is a 
stationary uniform product measure provided that 

h\u) (g) \u) = \u) (g) |a;) - \x) (g) 
Bi\u) — \x) + a\u), 

BL\u)^-\x)~a\u), (19) 

where a is a constant, which depends on the reaction rates. Defining the vector 
l^> by 

|7^) ~ h\u) (g (20) 
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then one obtains 



i?i — i?5 — Rg — 0, R2 + i?4 — 0, 

R^ + R^= 0, i?e + i?8 = 0, (21) 

and 



R2 i?B R?, 



U1U2 

where Ri is the ?th element of the vector \TZ). Defining p by 



(22) 



p:=a + b (23) 

(PI]) may be written as 

-w+i(l - p)2 + (cjia + wi4)(l - p)a + (cjia + ojnjil - p)b 
+ (wi6 + i^w)ab + cjisfl^ + cjigfo^ = 

^^51(1 - Pf + (t^52 + t^54)(l - P)a + (^53 + W57)(l - p)b 

+ (^56 + a;58)a& - ^+50^ + ^59^^ = 

W9l(l - p)^ + (^^92 + W94)(l - /9)a + (^93 + W97)(l - p)b 
+ (^96 + W98)a& + W95a^ - W+9&^ = 

{l021 + W4i)(l - p)^ + (-W+2 + t^24 + W42 - W+4)(l - P)a 
+ (W23 + W27 + W43 + '^47)(1 - P)b + (^25 + ^45)0^ 
+ (W26 + '^28 + ^^46 + ti^48)a^ + (^^29 + ^49)^^ = 

(a;3i + a;7i)(l - pf + (^32 + W34 + W72 + a;74)(l - p)a 
+ (-a;+3 + W37 + a;73 - a;+7)(l - p)b + (^35 + ^75)0^ 
+ (^36 + a;38 + a;76 + ^73)06 + (^39 + a;79)6^ = 

(cjgl + W8l)(l - pY + (t^62 + ^64 + W82 + a;84)(l - p)a 
+ (^63 + ^67 + ^83 + '^87)(1 ^ P)b + (^35 + ^85)0^ 

+ (-W+6 + ^68 + ^86 - t^+ii)ab + (^69 + ^89)6^ = (24) 

These constraints guarantee existence of a uniform product measure as a sta- 
tionary state. Depending on the model, the stationary state may be unique. 
For such models, for any initial state, this product measure is the final state. 



4 Product shock measures 

Let us assume the occupation probability of particle A (B) in the first m sites is 
oi (61) and in the latter (L — m) sites is 02 (62)- This is a single product shock 
measure, and the state vector is 



\u 
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where 




We briefly call it a shock measure. The shock state vectors \em), m = 0, 1, • • • ,L 
form L + 1 vectors and are a closed set under time evolution, if the evolution 
equations are in the form of 

H\e„i) = di\ern-i) + c?2|em+i) - {di + d2)|em), 1 < m < L - 1, (27) 

and 

H\eo) ^D[\ei)~D[\eo), 

H\eL)=D^\eL-i)~D'^\eL). (28) 

The shock position can perform one-site jumps with the rate c?i (^2) to the left 
(right) and at the left (right) boundary with the rate D'^ {D'2) to the right (left). 
As can be seen, the initial shock state can evolve into a linear combination of 
the shocks. It is said the shocks make an invariant sub-space. di,d2,D'i, and 
D2 are all non-negative and called the hopping rates of the shock position. 

A shock is immobile when all the hopping rates di , c?2 , D[ and Z?2 become 
zero. Then the shock does not move in the lattice, and the two stationary 
uniform product measures |eo), {cl) are stationary. 

4.1 Classification of shocks containing the occupation 
probabihties different from and 1 

When there is a shock in the lattice, it is divided into two parts. The occupation 
probabilities for particles A and B in each part are ak and bk , consequently < 
Cfc , bk < I- We also have < pfe := -|-5fc < 1. To have a shock, at least one of 
the particle densities in the first part of the lattice, oi or bi, should be different 
from that in the other part. Here we have assumed that ai ^ a2, 61 7^ 62, and 
Pi 7^ P2- A shock measure is constructed of two stationary uniform product 
measures, so the reaction rates should satisfy two sets of the equations (|21l22p . 
To classify the shocks, we should first define equivalent models. One may obtain 
models which apparently seem to be different, but there exist transformations 
which relate them together. We take these models as equivalent models. 

• Under the left and right exchange, some different models transform to 
each other so are equivalent. 

• The models which change to each other through the transformation AoB, 
are equivalent. Besides, because there is not any difference between the 
nature of an empty site and a site occupied by a particle, the empty site 
can be considered as a particle with the occupation probability (1 — p). 
Therefore, relabeling the single-site states does not lead to new models. 
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Considering the above mentioned transformations, classification can be done. 
Using the three parameters — , — , , and the two sets of equations 

0-2 h (1 - P2) 

(Plj) . two distinct models can be obtained. The two cases are 

. ai hi (1 -pi) 
02 02 (I-P2) 

^^)-<V<^T^, (29) 
a2 02 (I-P2) 

So it is enough to study the cases Ki, and K2. 
4.1.1 Ki 

Resulting from the condition — = -i , the non-zero reaction rates are 

02 02 

{a;23, 1^24, ^^27}, {^56,^58,^59}, 

{a;32, W34, W37}, {a;65,W68,'^69}, 
{a;42, W43, W47}, {w85,a;86, W89}, 

{ W72 , W73 ,^74}, { W95 , W96 , } , (30) 

in which particles A and B can convert to each other, but no particle can annihi- 
late to an empty state. Thus the total number of empty sites and consequently 
the total number of particles are conserved in the bulk. This model is symmetric 
under the the transformation A-H>B. As it will be seen the relations we obtain 
reflect this symmetry. We have 

(a;56 + a;58)ai6i - uj+r^al + aJsgbj = 0, 
(wge -I- a;98)ai6i + ^950? - w+96^ = 0, 



^32 + (^34 + ^^72 + ^^74 _ 
W23 + 1^43 + (^27 + "^47 ai 



(31) 



di is given by 



, ai 61 ai ai 

dl = (a;24 + W27 ) = (t^37 + W34 — ) = 

02 ai a2 oi 

= (^^42+^^43 — ) = -; (^^73 + ^^72 ) , (32) 

1 - P2 ai I - p2 bi 

dl and ^2 are related to each other through 

dl d2 



(1 - pi)ai (l-/02)a2' 
The boundary rates together with the reaction rates satisfy these relations 

/3l((w32 + W72)ai - (i:J23 + ^^43)^1) = 



(33) 



7 



= ("12 - aD^l^l - (-"32ai + "23^l)Pl' (34) 

di and ^2 can be described by boundary rates through 

-("21 + "3i)(l - Pi) + "12^1 + ai^bi 



di - 



Pi - P2 

("21 + "Sl)(l - P2) - "12^2 - "^362 



Pi - P2 

and also 

J/r, f/r 



(35) 



"21 bi — a^i ai, (36) 



Defining af^^ ■= afj + a^^, one arrives at 



^ i+r i+r\ I i+r . l+ri \ P2 /q^n 

("12 - "13 m = (-"aJ 02 + "23 02)7- (37) 



And finally the rates D'l , D'2 are 



J ^2 , "21 , Q2 -Qi , afzOi + 0^361 

— 0,2 - di- — \ — a2 \ r , 

h ai a2 (l-pi)pi 

D'^ = d,-d,^ + ^^d,^^^ + ^^^^^^±4^. (38) 

02 0,2 ai (1-P2)P2 



Setting all the transition rates equal to zero except for the diffusion rates 
1^737 '^37, ^^42, W24 and the exchange rates wsgjOJeg, the results of [13^ can be 
obtained 

^ ^2 J (1 ^ P2) 

^68=^^86, W24=W37 = ai , a;42 = W73 = di — (39) 

ai (1-Pi) 



(l-pi)ai (1 - P2)a2 



e/r, llr 
"21 "1 = "31 '^1' 

("12 - "13 )«iOi = (-"32 ai + "23 Oi)pi, 

-a|i(l -P i) - "3201 + ("13 + "23)^1 
«i = 



(40) 



61 -5: 



2 



, "31(1 - P2) + QS2Q2 - ("la + "23)^2 .... 
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, s Qf^i I \/-, N ai3&i + a^oai 

Z?! = (W42 - W24)P2 + ^ = (1^24 - ^A2){^ - P2 ) H 7^ ^ , 

61 (l-pi)pi 

- - C.42)P1 + ^ = (^42 - ^24)(1 " Pi) + "13^^ + (^^^ 

It is seen that the rates of diffusion for both particles A and B are the same as 
obtained in [TT] for a single-species model. So ASEP is a special example of the 
case K\. 

The shock will be immobile {d\ — di = D'^ ~ D'2 ~ 'ii) , \i 

W24 , "^42 , 1^27 J W72 , a;34 , a;43 , ^37 , 1^73 = (43) 
These reactions are 

0A o A0 
0A o B0 
0B o B0 

0B o A0. (44) 

The above transition rates are responsible for particle transportation from one 
site to another one. Whenever these rates are nonzero, the shock position moves. 
Thus, to fix the shock position, these rates should be zero. We also have 

"12 ' "13 ' "21 > "31 = (45) 

These boundary rates should be zero to maintain the total number of particles 
in the bulk. So we have 

W32 + ti^74 _ i>i 

W23 + ^47 ai 
W32ai - W23^i = a|2ai - "23^1 

("32 + a32)«2 = (a23 + "23)^2 (46) 
The other relations remain unchanged. 

4.1.2 K2 

In this model the occupation probabilities should satisfy the following relation 

(l-Pi)Qi _ (1 -P2)a2 , , 

hi - bl ^^'> 

Then the non-zero transition rates are 

W24, 1^29, 1^42, W49, CLI92, UJg4, (^37, 0^73, OJqs, LOsq (48) 

Therefore, besides the diffusion (cj24, W42 , W37, CJ73) and exchange processes 
(ojQs^ujse), provided that (j47]) is satisfied, the processes 0A o BB and A0 o BB 
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are also allowed. The occupation probabilities should also satisfy the following 
relation 

ai(l - P2) - a2(l - Pi) _ 0162 - a2^i . .„s 

bi{l - P2) - &2(1 - pi) ai(l - P2) - 02(1 - pi) 

This relation should be held in order to guarantee the possibility of the simul- 
taneous existence of the two uniform parts in the bulk. Then we have 

W37 = 1^86; 
'^73 = 

(w92 + a;94)(l - pi)ai - (^29 + ^49)^? = 0, (50) 



di and d2 are connected through 

di d2 



(1 - pi)ai (1 - p2)a2 



(51) 



di is given by 



ai (1-pi) ai bi ai-a2, 

"1 = ^^37 — Ti \'^73 — 1^24 +t^29 ) = 

0-2 (1 - P2) a2 ai Pi - P2 

{I- Pi), , b2ai-a2, (l-pi)ai a2 «lYro^ 
— -(W42 + W49 ) = 7 — r rT('^92T a;94— {52) 



(1 - P2) 02 pi - P2 (ai&2 - a2bi) 62 bi 

Defining af, + a^, =: a^.f'', then D[ , DL will be 

^'1 = 6^(-"3r (1 - P2) - a^3r«2 + (air + "2r)&2) 
^2 = ^-^(a^r (1 - pi) + «3rai - («!r + «2r)''i) (53) 

O2 — Oi 

The boundary relations have been included in Appendix. A question which 
may arise is that: is there any solution for these set of conditions? These con- 
ditions are some constraints on the reaction rates, and occupation probabilities 
of the particles, and 6^, which should be checked. In addition to these equa- 
tions, there arc also some inequalities; reaction rates should be nonnegative and 
< ai,6i < 1. It is difficult to check all these analytically. We have solved 
numerically these set of equations and inequalities. It is seen that there are 
solutions for these set of equations. So all the conditions on the parameters can 
be fulfilled simultaneously. 



This case does not have a similar one-species model. To have a one-species 
analog, the processes in which all the three single-site states are involved should 
vanish, i.e. these processes A0 o BB should not happen. This means that 
just the diffusion and exchange transition rates can be held which leads to a 
contradiction to the inequality — < ^ < tt— 

Moreover, immobile shock is not possible in this case. If we set all the 
hopping rates of the shock position equal to zero, then all the transition rates 
become accordingly zero. So, this model does not admit immobile shocks. 
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4.1.3 Generalization to more than two species 

Let us consider an 71-species system. The local hamiltonian h, acting on two 
adjacent sites, is an (n + 1) x (n + 1) matrix. We assume a single-shock measure 
consisting of the single-site state vectors \u) and \v). It can be seen from the 
equations PHIET]) . if the equations governing the evolution oi h\u)(^\u), /i|u) (g) 
\v) and h\v) ^\v) are given, the dynamics of the shock in the bulk of the lattice 
is fully determined. To facilitate the calculation, we choose a new basis like 
C = {\u), \v), \xi), . . . , |x„_i)}, where \xi)^s are arbitrary vectors ineffective in 
the final results. Adopting the same strategy as Section 3, from expanding the 
three vectors h\u) \u),h\u) ^ \v) and h\v) ^ \v) in the new basis and then 
substituting these expansions into (j27p . one obtains the following equations for 
the uniform parts of the shock 

n-1 

h\u) (g) \u) = Ei{\u) \v) - \v) (g) \u)) + Ef,+i{\u) d) \Xf,) - \Xf,) «) |u)), 

M=l 

n-1 

h\v) (g) \v) = -Fi{\v) (g) \u) - \u) (g> \v)) + J2 \x^,) - g) \v)), (54) 

and the following equation for the shock position 

h\u) (g \v) = (rf2 - Fi){\u) g) \u) - \u) (g \v)) + {di - Ei){\v) (g \v) - \u) (g \v)) 

n — 1 n — 1 

fi—l M— 1 

E's and F's are some constants depending on the reaction rates and probabili- 
ties. When all the occupation probabilities are different from zero and one, the 
equations (|54p seem to be sufRcient to determine the non-zero rates. In this 
case, the equation (I55t for the shock position either manages to link the two 
uniform parts or completely fails and rejects the validity of the whole model. 
The equations ([M)) are both in the following form 

h\w) g) \w) = \w) g) \x) - \x) g) \w) (56) 

We define \Ti) := h\w) (g) \w) whose elements are i?fc's. We have obtained some 
relations for these elements by generalizing our previous knowledge of what 
these relations are for one and two-species systems, (|2T|) . We have 

i?(^_i)(„+i)+^ ,^^l,...,n + l (57) 

These relations belong to the two-site states in which either both sites are empty 
or occupied by particles of the same type, for example 00 or AA . We also have 

-R(V-l)(n+l)+0 + -R(0-l)(n+l)+i/. = 0' S ^ = 1, ■ ■ • , (58) 

[ = 2,...,n + l 
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Each relation of this kind is related to those two-site states which are the 
mirror of each other, for example A and A 0. 

It is known that, in the case of single-species models there is only one model 
with < p < 1, admitting shock as a random walker [TT]. This model is the 
asymmetric simple exclusion process (ASEP). For the two-species case, we have 
found two models here. One of these models is a generalization of ASEP, and 
the other one has no single-species analog. 

Let us consider a three-species model. The single-site states respectively are 
empty site denoted by and occupied sites by A, B and C. Then the two-site 
states become: 



1) 00 


5) A0 


9) B0 


13) 


C0 


2) 0A 


6) AA 


10) BA 


14) 


CA 


3) 0B 


7) AB 


11) BB 


15) 


CB 


4) 0C 


8) AC 


12) BC 


16) 


cc 



and the equations ([57]) and (1551) give 



i?2 + ^5 — Oi ^3 + ^9 — Oj ^4 + ^13 — 

Rt + Rio =0, Rs + RiA - 0, i?i2 + i?i5 - (60) 

To do classification, one should compare the following quantities with each other 

{I -Pi? of 6? c? (I-Pi)ai 
{I-P2Y ' ai ' bl' 4' (1-P2)a2 ' 



(1 - pi)bi (1 - pi)ci aibi oici bici 

(1 - ^2)^2 ' (1 - P2)C2 ' 02^2 ' a2C2 ' ^2C2 

where a, b and c are the occupation probabilities of particles and p is defined 
as {p := a + b + c). These parameters are suitable criteria for the ratio of 
the number of the two-site states of the two uniform parts. As we know, how 
much greater the number of a definite two-site state is, the associated process 
is more likely to happen. Therefore, the reason for which the parameters (|6ip 
are helpful to introduce tractable models is that they can be used to control the 
processes in the lattice. Since the parameters (j6ip are the square or product of 
the parameters i 7 and , it will be easier to first compare these 

parameters. 

Three distinct models will be obtained for a three species lattice. 

• The cases with one inequality like 

ai bi ci (1 -Pi) , s 

0-2 02 C2 (I-/02) 
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For this case, the non-zero transition rates are 



i,j = 6,7,8, 10,11,12,14,15,16, 
ujik, l^k, Z,fc = 2,3,4,5,9, 13. (63) 

in which particles A, B and C can convert to each other, but the number 
of empty sites is conserved. Resuhing from (I62p . this model is invariant 
under the transformations A B, A C and C O B. It shows that, 
from the four single-site states of this case, the states A, B and C are 
similar to each other. Another model of this type is 

- = ^ < - = (64) 

0-2 02 C2 (I-P2) 

The non-zero transition rates are 

We/, e^/, e,/ = 2,3,5,8,9,12,14,15, 
^gh, 9 7^h, g,/i = 6,7, 10,11, 

o^p, o,p= 1,4, 13, 16. (65) 



This model has the symmetries A -o- B and C -O- 0, so the states A with 
B and empty site with C are similar to each other. Therefore, for the 
cases with one inequality, one expects an analogous single-species model, 
resembling two reduction in the number of particles. 

• The cases with two inequalities like 

ai h ci (1 -pi) , , 

— = — < — < r (66) 

0-2 02 C2 (I-/O2) 

This model has the symmetry A o B, thus it is expected to have a similar 
two-species model, resembling one reduction in the number of particles. 
The nonzero transition rates are 

W4,13 , 1^13,4 

w.'j', iV/, = 6,7,10,11, 

uji,^k', I'^k', r. A:' = 8, 12, 14, 15, 

e'y^/', e',/' = 2,3,5,9,16. (67) 



e'f, 

in which, besides the diffusion and exchange processes, particles A and B 
can convert to each other and the processes (0A, A0, 0B, B0 ^ CC) are 
also possible provided that the following relation is satisfied 

^^-^ = 4 (68) 
(1 - P2)a2 ci 

There is another essential condition for the occupation probabilities which 
is 

((ai 5i)(l - P2) - {02 + h2){\ - pi)f = 
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= (ci(l - P2) - C2(l - /9i))((ai + 6i)c2 - (02 + ^'2)ci) (69) 
This relation comes from the equation (|55|) for the shock position. 
The cases with three inequahties hke 

< ^ < < (70) 

are expected to be originally three-species models. The nonzero transition 
rates are 

'-^7,10 , "^10, 7 J '^4,13 , "^13,4 , 

t^a'h', g'y^h', 5'>' = 2,5,12,15, 
Wo'p', oVp', o',p' = 8,11,14, 

g'^s', g',s' = 3,9,16. (71) 

with the conditions 

(l-pi)ai bici (l-pi)&i cl flici 6^ 



(1 - P2)a2 ^2C2' (1-^2)^2 C2' a2C2 6: 



(72) 



In this case, besides the diffusion and exchange processes, resulting from 
the above relations, the processes (0A, A0 o BC, CB), (0B, B0 o CC) 
and (AC, AC -H> BB) are also allowed. The occupation probabilities 
should also satisfy another relation coming from the equation (|55p. 



Accordingly for n-species systems, it is expected to find n distinct models. One 
for which there is a similar single-species model, one with a two-species analog, 
• • • , and finally there is an n-species model which has no fc-species analog with 
k < n. 



4.2 shocks containing the occupation probabihties and 
1 

Now, we study the shocks in which some parts of the lattice is empty (or equiv- 
alently completely occupied). 



4.2.1 pi^0,l,p2 = 

Assume 02 = 62 = and oi, 61 7^ 0, 1. In this case, loh = to prevent particle 
production in the empty part of the shock, and subsequently lou = 0, otherwise 
all the two-site states can convert to 00 but there is no way to leave this state. 
All the other remaining rates can be non-zero. Besides the relations ([21]) for 
pi,bi,ai, the other relations are as follows. The parameters di, ^2 are given by 

di = —[{uj42ai + i^4:3bi)ii ~ Pi) + {(^46 + ^is)aibi + u^zaf + uj4gbl] 
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-^[(w72ai + a;736i)(l - pi) + (wre + a;78)ai6i + ivr^al + W79?>i] 
-(w24ai + W2761) = 7j — , (a;34ai + W3761) 



(l-pi)ai (l-pi)6i 

= -o (W54ai + W5761 ) = — (W64CH + W67&1 ) 

d]^ dlOl 

= -r^{u)84ai + ivsrh) = ^{uigAai + ujgrbi). (73) 

Oifli 

The boundary rates together with the reaction rates satisfy these relations 
(ai2 - airi)aibi + (03201 - Q23''i)Pi = (-w+4ai + W4761 + d2ai)pi = 

= [{UJ62 + W64)ai + (t^63 + W67)&l](l - Pi) + ^650? + W69&1 

+(-a;+6 + W68)ai^i, (74) 
and we have 

- d2Pi(l - pi) + dipi = -(a2i + a3i)(l - Pi) + Q!i2«i + oei^h, (75) 

Particles should not enter the right boundary, so 0:31 = = 0. The relations 
between the boundary rates are 

02161=03101. (76) 
Using the definition a^j~'^ := afj + a^j, one arrives at 



(air - ai+naibi = (-a^g^ai + ^^^Opi- (77) 



D[, D'2 are 



D; = Ji-d.(l-pi)+ ''"'''+°"'" . (78) 

pi 

The shock will be immobile if all the transition rates from the fourth and seventh 
rows and columns of h except for become zero, i.e. the processes starting 

from (ending to) the states A0, B0 which make the shock position move to the 
right (left). Also the following rates at the boundaries should be zero 

"21 = "31 = "l2 = o^ia = ot\2 = (Aa = 0- (79) 

which means no particle enter or leave the boundaries. The other relations will 
remain unchanged. 

4.2.2 pi = l,p2 = 

Consider 02 = 62 = and oi + 61 = 1. We assume oi, 61 ^ 0, 1 to have a two 
species model. The non-zero transition rates are 
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Wl2 


Wl3 
W23 


UJ14 






Wl7 






W32 
















UJ42 


^^43 








lx)47 






OJ52 


W53 


W54 






W57 


W58 


W59 


^62 


^^63 


W64 


^65 




a;67 






LUr2 


W73 


W74 












W82 


^^83 


W84 


W85 


W86 


W87 




^89 


W92 


^93 


<^94 


W95 


t^96 


W97 


^98 





There are no constraints on uji2 and uj^. One should expect it, because the 
states (0A) and (0B) do not exist in the final state. We also have 



(W56 + W58)ai&i - W+sOj + W5961 = 

(a;96 + a;98)ai&i + ^95^1 - w+961 = 0. (80) 



di , d2 are given by 



di = Wi4ai + ^1761 

d2 = -2 (W54ai + ^bjbl) —(^6401 + t^67&l) 

af fliOi 

= 7 (W84ai +W87&1) = 70(^94(11 +W97^l)- (81) 

The following relations should also be satisfied by reaction rates 

a32ai - a23^i = (-'^+6 + W68)ai6i + wesa? + t-Jeabf 

= ((^2 + rfi)ai - a;+4ai + W4761. (82) 

The left part of the shock should be fully occupied, so a{2 = = and the 
right part of the shock should be empty, thus a2i = a^i = 0. The relations 
between the boundary rates are 

a^2ibi = Q-iiai 

ial2 - Q!i3)ai&i = -("32 + Q!32)ai + ("23 + "23)^1 (83) 
D[ , D2 are given by 

D2 = al2ai+alsbi. (84) 

The shock will be immobile if we eliminate the processes (A0, B0 — >• 00), i.e. 
Wi4 = CJ17 = 0, to prevent the shock position from moving to the left, and also 
the processes (A0, B0 — AB, BA, AA, BB) should vanish to prevent the shock 
position from moving to the right, i.e the transition rates W47,a;74 are the only 
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non-zero rates of the forth and seventh columns of h. The fohowing boundary 
rates should be also zero 

"21 = "31 = "^2 = "^3 = 0. (85) 

which means no particle enter or leave the boundaries. 

Some other shapes of the shocks for two species systems can be also predicted 

(ai=0, 02^0,1, 6i^0,l, 62 = 0) 

(ai=0, a2^0,l, h = l, 62 = 0) 

(ai^O,l, 02^0,1, 61 0,1, 62 = 0,01+61^0,1) 

(ai7^0,l, 02 7^0,1, 617^0,1, 62 = 0,01+61 = 1) (86) 

It can be shown that the models with oi =02, (oi, 02 7^ 0, 1), 61 7^ 62, (61, 62 7^ 
0, 1), and pi, p2 7^ 0, 1, do not exist. 



Dynamical phase transition for two-species mod- 
els 



In this section dynamical phase transition for the three cases Ki, (pi 7^ 0, 1, p2 = 
0), and (pi = 1,P2 = 0) is studied. In [17 , phase transition in single-species 
models possessing shock solutions is studied. It will be shown that there are 
three phases. For some region of the parameter space, the relaxation time is 
independent of the reaction rates at the boundaries. Changing continuously 
the reaction rates at the boundaries, there is a point where the relaxation time 
begins changing, so that at this point there is a jump in the derivative of the 
relaxation time with respect to the reaction rates at boundaries. This is the 
dynamical phase transition. So the dynamical phase transition studied here is a 
discontinuity in the derivative of the relaxation time (from zero to nonzero) with 
respect to reaction rates in the bulk and at the boundaries. We use the method 
presented there, for our two-species models. Defining the two parameters A, 
and B as 



A 



di - D' 



\/did2 



B 



d2-D[ 



(87) 



it is shown in [T7] depending on the phase of the system the relaxation times of 
the models may be Tq, T^h or Tg 



To 



di - D'. 



{Vdi - Vd2y 

Three phases may occur 



D'2{di 



D'2 



1) A<l,B<l 

2) A<l,B>l 

3) A>l,B<l 



Tk 



To, 
Ta- 



D[id2-D[-di)- 



(89) 
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The case ^ > 1 , > 1 do not occur in none of the above mentioned models. 
Let us first consider the case Ki. we have 



^ (W24 + W27^)(W42 +^^43^) (^24 + W27^)(W42 + '^43^^) 

^_ (0^24 +0^27^)-^ _ (w42+W43^)-^^q;7)^ 



/(W24 +a;27|7)(W42 + W43I7) ^ (W24 +^27^) (^42 +^43^) 

If (w24 + W27^) > (w42 + ^43^) One arrives at 



w (^24 + W27 — )(W42 + W43 — ) > (a;42 + W43 — )• (90) 

V fli di di 

So ^ (uj2i + W27 ^ ) (W42 > (W42 + W43|^) - which gives A < 1. For 

(W24 +W27 — ) > (W42 + W43— ), (91) 

cii ai 
it can be shown that B is also less than one. If 

(W24 + W27 — ) < (W42 + W43 — ) , (92) 

fli tti 

one again can show that both A, and B are less than one. Thus this model 
belongs to the region 1. So the model Ki possesses no dynamical phase tran- 
sition, and the rcilaxation time is Tq. To has no dependence on boundary rates 
and only depends on the reaction rates in the bulk. This result is the same as 
the one found for ASEP, which is a single-species model. It should be noted 
that this model is a generahzation of ASEP. If one forgets about the difference 
of two particles, Ki changes to ASEP. 

In the case of pi 7^ 0, 1, p2 = 0, defining 

fll := — 7K'^42 + W43 — )(1 - pi) + (W46 + W48)&1 + W45ai + W49— 1, 

(1-pi)'-' oi ai^ 

0.2 ■= (W24 + W27 — ), 

ai 



one arrives at 



a[2(ii + (y[J>i 



VL2 

A= ^ 



g ^ - Pi) ^ i'^-Pl) f>l /ggN 
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if fii < ri2, then B < 1 and A can be smaller or greater than one and if fii > 
A < 1 and B can be smaller or greater than one. So, in this case the regions 
1,2 or 3 are possible. There are three distinct relaxation times. Tq is bulk 
dominated and depends only on reaction rates in the bulk, while T4 and Tg 
depend both on the reaction rates in the bulk, and at the boundaries. So in this 
model there are three distinct phases available for the system, and this model 
may experience dynamical phase transitions. Changing the reaction rates at the 
bulk or at the boundaries may lead to phase transitions. 
And finally, in the case pi = 1, p2 — 0, defining 

n[ := (i:Ji4ai + ujijbi), 

^2 ("^54 + 1^64 + + a;94)ai + (1LJ57 + ujej + ujg,7 + ^^97)^1 , (94) 

it can be shown that this model's parameters may be in any of the regions 1, 
2 or 3. So this model may also experience dynamical phase transitions. It is 
known that single-species models have three distinct phases [T7]. Here, It is seen 
that the two-species models considered here are also restricted to the regions 
1,2, and 3. 

6 Discussion 

In this work, we presented a method for classifying n-species particle systems, 
possessing single shock solutions. This classification provides a simple strategy, 
adopting only the occupation probabilities of particles. We applied this method 
to two- and then three-species systems. 

We have made a detailed study on the four two-species models Ki, K2, 
(pi 7^ 0,1, p2 = 0) and (pi = l,p2 = 0). For all these models some new 
results have been found here. In [16] the three models Ki, [pi ^ 0,1, y02 = 
0) and {pi = l,p2 = 0) have been studied. In this article, with the same 
assumptions, more general solutions have been found. This shows there are 
some unnecessary constraints in |16| . which could be eliminated. For example, 
for (pi = l,/32 = 0), there should be no constraints on uji2, uji^, because the 
states 0A, 0i? do not exist in the final state, and also the transition rates 
1^561 W58, W59, wgS) t^eS) ^69, L0s5, Wge, wgg, W95, ^96, 1^98 from the 5th, 
6th, 8th, 9th columns of the local hamiltonian h were set equal to zero, while 
generally they could be nonzero (see subsection 4.2.2). 

The model investigated in [TS] resembles the model K2 of this paper in some 
of the initial assumptions; the nonzero bulk transition rates and the condition 
(|47p of this paper are the same. However, each paper has been presented a 
completely distinct solution. In [15] . the occupation probabilities have been 
assumed to be equal to and 1 in one part of the lattice, which required three 
of the ten transition rates to be zero. This model is a particular example of 
the more general model pi ^ 0,l,p2 = (Subsection 4.2.1). We have found 
another solution with the assumption that the occupation probabilities should 
be necessarily different from and 1. Consequently, one should note that the 
two papers presented two different models. 
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Moreover, the model studied in Subsection 4.1.1 with a special choice of 
symmetry has been previously studied in [M]. A special example of this model 
when the reaction rates are eliminated has also been presented in [T3]. We de- 
rived the same relations for this case in Subsection 4.1.1. 

Acknowledgement: A. A. was partially supported by the research council 
of the Alzahra University. 



7 Appendix 

The simplest forms of the boundary relations of the case K2 are 



1) 




- P2) 


+ ("12 


+ 0^2)02 


- (223^2 + 










- P2) - 


1- a^^^a2 


- ("13 + 


"23)^2) , 
O2 


- ai 
-61 


= /i 


2) 




- P2) 


+ {a[2 


+ "32)«2 


- "23^2 + 










- P2) - 




- ("^3 + 


r \u \ 
"23)02) 


- ai 
-bi 


= -Ii 


3) 


"«2l(l 


-Pi) 


+ ("L 


+ "32)«1 


- "23^1 + 










- Pi) - 




/ i , I \i \^2 
("13 + "23)^1) ^2 


- ai 
-61 


= h 


4) 


-«2l(l 


~Pl) 


+ {a[2 


+ "52)ai 


- "23^1 + 










-Pl)- 




- ("I3 + 


r \i \^2 
"23)01) . 

O2 


- fli 
-bi 


= -l2 



Defining 



5) 7— -T (-"31(1 - P2) - "32^2 + (0^3 + 0^23)^2) = /3 
02 — Ol 

6) 7-(a^3i(l-Pi) + "32ai-("i3 + "23)&i) = ^4 (95) 

02 — Ol 



/ («l&2-fl2&l)(pi-p2) \ . . 

(ai-a2)(6i-52) ' ^^^^ 



-^1, ^3 1 -^2, ^4 become 

(ai&2 - a2&i)^ (/Ol - P2) 



/2 -dr 



010261 (fli - a2)(6i - 62) 



/4 = rfi(l + M^ 

ai(l - pi) 



/3 = di( I +M^ii-^). (97) 
b{ ai(l-pi) 



As it can be seen, if one of the constants di, /i, /3, I2, I a is given, the other 
constants will be obtained from the above mentioned relations. 
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